QCD prediction for the non-DD annihilation decay of -0(3770) 
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To clarify the marked difference between BES and CLEO measurements on the non-DD decays 
of the 0(3770), a l 3 Di-dominated charmonium, we calculate the annihilation decay of 0(3770) in 
NRQCD. By introducing the color-octet contributions, the results are free from infrared divergences. 
The color-octet matrix elements are estimated by solving the evolution equations. The S-D mixing 
effect is found to be very small. With m c — 1.5±0.1GeV our result is r(0(377O) — > light hadrons) = 
467^33g KeV. For m c = 1.4 GeV, together with the observed hadronic transitions and El transitions, 
the non-DD decay branching ratio of 0(3770) could reach about 5%. Our results do not favor either 
of the results of BES and CLEO collaborations, and further experimental tests are urged. 
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Heavy quarkonia decays play an important role in un- 
derstanding quantum chromodynamics (QCD)[l|]. These 
include not only the determination of the running strong 
coupling constant a s from S-wave decays J/ip —> ggg and 
T — ► ggg, but also the study of factorization from the 
P-wave annihilation decays, where appear infrared (IR) 
divergences in 1 P\ — * ggg and 3 Pj — ► gqq @, A tra- 
ditional way to treat the IR divergences was to use the 
quark binding energy or the gluon momentum as cut- 
off to estimate these IR divergences, but this is model 
dependent and breaks factorization of short and long 
distance processes. In a new factorization scheme 
was proposed to absorb the IR logarithms by new non- 
perturbative parameters, the color octet matrix elements. 
Based on the non-relativistic nature of heavy quarkonia, 
an effective theory, Non-Relativistic QCD (NRQCD) was 
developed Q, in which the inclusive annihilation decays 
can be calculated in a systematic way by double expan- 
sions in terms of a s and Wjthe relative velocity of quarks 
in heavy quarkonium. In [(J, 0, S| , the authors calculated 
QCD radiative corrections to the light hadron (LH) de- 
cays of P-wave charmonium in NRQCD, and showed ex- 
plicitly the cancelation of infrared divergences at the next 
to leading order (NLO). In j| a more complete and pre- 
cise NLO calculation for the P-wave decay perturbative 
coefficients in NRQCD is given (see also [10(| ) . At NLO in 
a s , the NRQCD predictions for the relative decay rates 
°f XcJ LH are consistent with more updated data (see 
Chapter 4 of (H). Moreover, the relativistic corrections of 
S and P-wave electro mag netic quarkonium decays have 
been given at order -l? 7 [llj] ■ As for the D-wave, in 12, 13] 
calculations of 3 D,/ — > ggg decays were given but suf- 
fered from IR divergences; while in [§J only the leading 
order (LO) color-octet contribution to 3 D,/ — ► LH was 
given. So for the D-wave a complete calculation for the 
IR cancelation and radiative correction in NRQCD is ap- 
parently needed. 

Phcnomcnologically, for the 3 Di (J PC = 1") char- 
monium state -0(3770), there is a long-standing puz- 
zle in its non-DD decays that the 0(3770) might have 
substantial decays not into D°D° and D+D~. BES 



earlier reported two results based on different analy- 
sis methods: Br(0(377O) -> non-DD) = (14.5 ± 1.7 ± 
5.8)%fll, and Br(0(377O) -> non-DD) = (16.4 ± 7.3 ± 
4.2)%[l|. In contrast, CLEO[l3 measured the cross 
section a{e r e~ -> 0(3770) -> non-DD) = -0.01 ± 
0.08^Q'gJnb. Very recently, with the first direct mea- 
surement on the non-DD decay, BES gives a(e + e~ — > 
0(3770) -► non-DD)_ = (0.95 ± 0.35 ± 0.29) nb and 
Br(0(377O) -► non-DD) = (13.4 ± 5.0 ± 3.6)%[l7|. Evi- 
dently, the two collaborations give very different results 
of the non-DD decay of 0(3770). Meanwhile, a number 
of experiments to search for the exclusive hadronic non- 
DD decays of 0(3770) have been done by BES [Hi and 
CLEO[l9j, but no significant signals are found. 

At least two kinds of non-DD decays of -0(3770) have 
been observed. The hadronic transitions 0(3770) — > 
7r + 7r~J/0 was first observed by BES with a branching 
ratio of (0.34±0.14±0.09)%[2l). and was later confirmed 
by CLEO with a somewhat smaller branching ratio 
Br(0(377O) -> tt+tt" J/ip) = (0.189±0.020±0.020)%[22j, 
and the ir°ir°J/iJj and rjJ/ip modes were also seen with 
each having a branching ratio of about one half of that 
of 7r + 7r~J/ip (22)]. These results are within the range 
of theoretical predictions based on the QCD multipole 
expansion for hadronic transitions [2JJL With the total 
width of 23.0 ± 2.7 Mev for V(3770)|2|, the width of 
all hadronic transitions is about 100— 150 KeV. Another 
kind of non-DD decays of ^(3770) are the El transi- 
tions 0(3770) — > 7 + Xcj (J=0,l,2), and their widths are 
measured by CLEO to be 172 ± 30, 70 ± 17, < 21 KeV 
for J=0,l,2 respectively [2ij], which are in good agree- 
ment with predicted values 199, 72, 3.0 KeV in a QCD- 
inspired potential model calculation with relativistic cor- 
rections [2H (see also (2fj|,[27j]). The width of all El transi- 
tions ^(3770) -► 7 + Xcj (J=0,l,2) is about 250±50KeV. 
The above mentioned hadronic and El transitions only 
contribute 350-400 KeV and 1.5-1.8% to the non-DD de- 
cay width and branching ratio of ^(3770). 

To clarify the puzzle of -0(3770) non-DD decay, in this 
letter we will give a complete infrared safe NLO QCD 
corrections to the annihilation decay rate of the 0^(3770) 
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in the framework of NRQCD. Since v 2 ~ a s (m c ) w 0.3 
in charmonium, the relativistic corrections are also im- 
portant and should be considered in the future work. 

The -0(3770) can be viewed as a 1 3 D 1 dominated state 
with a small admixture of 2 Si, and expressed as (see 

e.g. mm) 



1^(3770)) = cos 6>|l 3 L>i) + sin 6»|2 3 S*i), 



|0(3686)) 



sinflll 3 !?!) +cos6»|2 3 S* 1 ), 
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where 8 is the S-D mixing angle and it is about (12 ± 
2)° by fitting the leptonic decay widths of 0(3770) and 
■0(3686). Then the LH decay width of 0(3770) is 



T(0(377O) -> LH) = cos 2 6T(l 3 L>i -► LH) 
sin 2 6T(2 3 5i -► LH) + IF, 



(2) 



where IF stands for the S-D interference term. The cal- 
culation of S-wave decay at order a 3 and leading order 



in v 2 is trivial, and it gives 

r(2'^M)J^' a ^: 9 l (3) 



-in 



81m 2 



where i? 2 s(0) is the 2 3 S"i wave function at the origin. 
The S-D interference term IF in Eq.(2) is infrared finite 
at leading order in v 2 and a s , and can be obtained by 
combining the \ 3 D\ — > 3g with 2 3 5i — > 3g amplitudes 



IF = 2sin6>cos6>- 




324m4 



^(0),(4) 



where Rp{0) is the second derivative of the l 3 Di wave 
function at the origin. 

We now proceed with the calculation of the main part, 
the D-wave quarkonium LH decay. In NRQCD, the in- 
clusive annihilation decay of 3 D\ at leading order in v 2 
is factorized as 



T( 3 Dj -> LH) = 2Imf( 3 D l j ] )H m + ]T 21mf( 3 Pf ] )H P8J + 2Im/( 3 sf )F S8 + 2lmf( 3 S[ 1] )H sl , 

I 



(5) 
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where Im/(n) is the imaginary part of the QQ — > QQ 
scattering amplitude, and can be calculated perturba- 
tively . And the corresponding non-perturbative matrix 
elements are 

(g|0 1 ( 3 D 1 )|g) _ (H\Q 8 ( 3 Pj)\H) 

tlD\ = a ,WP8J — , 



(H|0 8 ( 3 gi)|g) „ (HIO^S^H) 
Hss = 5 >-"Sl = 5 > ("J 



I'll 7- 



where H is ip(l 3 Di). Those four-fermion operators of S- 
wave and P-wave are defined in Q , and here we only give 
the definition of the D-wave four-fermion operator (the 
normalization of the color singlet four-fermion operators 
agree with those in [9(): 



Oi( 3 -Di) = 



10iV r 



(T) 



where T l = S lj and S kl = {^) 2 {D 1 -\D 2 S ij ). 

We calculate the short distance coefficients at order a 3 , 
and details of our calculation will be given elsewhere. The 
S-wave and P-wave short-distance coefficients have been 
calculated in [§], and our calculated results agree with 
theirs. The D-wave short distance coefficients presented 
here are new, and they are 



2Im/( 3 Si 



3 g [ih_ 40a 3 (^ 2 -9) 
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2Im/( 3 sf) = ^(36N f n + a s (5(-657 + 67^ 2 ) 
lUo 

+ N f (642 - 207V/ - 27tt 2 + 72ln2) + lUfoNfh 
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„ rsl 5a 3 (4107- 647V f -4147r 2 + 487V f Zn^) 
2Im/( 3 pj 8] ) = — ^~ 
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774tt 2 + 1008/^2 + 1728/3 ln-£— + 240iV/Zn— ^)), 

(8e) 
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2/m/: ! /), : : = ' 
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where (3q = 11JVc 6 2 / , N c = 3, Nf is the number of 
flavors of light quarks, fi and /ia are renormalization 
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and factorization scales respectively. We consider ten 
processes to get the short distance coefficients in Eq.[8], 
including gg, ggg, qq, and qqg final states. The contri- 
butions of qq and qqg processes are labeled by the powers 
of Nf. 

After calculating the short distance coefficients, we 
come to determine the long-distance matrix elements. 
In the P-wave charmonium decay, at leading order 
in v 2 there are two four-fermion operators HI and 
#8(3, while in the case of D-wave, there are four 
independent matrix elements under the heavy-quark 
spin-symmetry. They are Hm, Hps, Hss, H$i, where 

tt 1 (H\Q 8 \ 3 P )\H) '_ 4(H\Q 8 ( 3 P^H) _ 20(H\O 8 ( 3 P 2 )\H) 
^PS ~ ^ 3^3 - ^3 , 

and these relations can be derived by considering the El 
transition from 3 Di to 3 Pj. In NRQCD, Hm is related 
to the wave function's second derivative at the origin, 
while for the other three, in the absence of lattice sim- 
ulations and phcnomcnological inputs, we will resort to 
the operator evolution equation method suggested in [B[ , 
where the authors give the result of the matrix elements 
in the P-wave decay. Here we derive the following matrix 
elements in the D-wave case 



Hps — — 



5 8C F a s (/jA ) 
9 3/3 a s (/ZA) 



Hs8 = 



CpBp 8 ^ 2 ^2/ gf(MAo) 

2 3/3 a s (p\) 



Dl, 



H, 



si 



Cp 8 

47V 3/V 



2| n 2j« S k) 



)Hm, 



(9a) 



(9b) 



(9c) 



TABLE I: Subprocess decay rates of D\ charmonium, where 
v 2 = 0.3, = 2m c ,a s (2m c ) = 0.249. 
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FIG. 1: Renormalization scale ^-dependence of the decay 
width of charmonium Z D\ to light hadrons. Here NLO means 
LO contribution+NLO correction 



where Cp = §, -Bp — ^ . We choose the region of validity 
of the evolution equation: the lower limit fj,\ = m c v and 
the upper limit ^a of order m c . 

With both the obtained short distance coefficients and 
long distance matrix elements, we predict the LH decay 
width of D\. The renormalization proceeds by using 
the MS scheme for the coupling constant a s and the 
on shell scheme for the charm quark mass. For conve- 
nience, we take the factorization scale /ia to be the same 
as the renormalization scale /i of order m c . We choose 
the pole mass m c = 1.5GeV, v 2 — 0.3, //a — m cV,^A = 
2m c , a s (2m c ) = 0.249, Nf = 3, A QCD = 390MeV, H m = 

15 'gffl 2 = 0.786 x lQ- 3 GeV^. At 0{a 2 ), the LH de- 
cay involves three subprocesses ( 3 Po)s - * gg, ( 3 ^2)s 
gg, ( 3 <Si)8 —> gg, and the decay width is estimated to be 



r( 3 £>i 



LH) = 0.205McV. 



(10) 



At 0(a 3 ), there will be seven more subprocesses 
( 3 Si)i, 8 -» ggg,( 3 Pi)s - ggg, ( 3 Pj)s - <?95,( 3 #i)i -» 
ggg involved, and the result turns to be 



Our result shows that in NRQCD factorization the NLO 
QCD correction is even larger than the LO result. The 
numerical values for all subprocesses are listed in Table 
1. If we choose fi\ = Tn c , a s (m c ) = 0.369, the values 
of LO and NLO (the sum of LO contribution plus NLO 
correction) become 0.28 and 0.68 MeV respectively. The 
renormalization scale /i dependence of the decay rate is 
shown in Fig.(l). We see that the /z-dependence at 0(a 3 ) 
is rather mild when /i > 0.9m c . For simplicity we take 
(i = 2m c , where the logarithm term In^r- = 0. 

With the pole mass m c = 1.5GeV,cv s (2m c ) = 0.249, 
\ r id(°)\ 2 = 0.015GeV 7 , and |i? 2 s(0)| 2 = 0.529GcV 3 , 
9 = 12°, we find that the three terms on the right hand 
side of Eq.(2) contribute 417, 5.3, 44 KeV respectively to 
the LH decay of "0(3770), and result in 



r(V>(3770) -> LH) = 467KeV. 



(12) 



r( 3 £>i -> LH) = 0.436MeV. 



(11) 



Our result shows that the D-wave LH decay is dominant, 
and the S-D mixing only has a very small effect on the 
-0(3770) LH decay. One important uncertainty of our 
prediction is associated with the long-distance matrix el- 
ements, especially the color-octet matrix elements. Using 
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the same evolution equation method in \c.j decays, we 
find the ratio of color octet to color singlet P wave de- 
cay matrix elements agree with the lattice calculation [2^| 
to within about 20% and with the phenomenological 
values 0, G3 to within about 30%. This might indicate, 
though not compcllingly, that the uncertainty related to 
the matrix elements calculated using the evolution equa- 
tion in the D-wave decays are also about (20-30)% or 
(with more confidence) less than 50%. Other uncertain- 
ties such as the relativistic corrections and higher order 
QCD radiative corrections are beyond the scope of the 
present study. On the other hand, however, we find the 
decay rate to be sensitive to the value of the charm quark 
mass. If we choose the pole mass m c = 1.5 ± O.lGeV, 
a s {n) — a s (2m c ), and fix other parameters as before, 
then our prediction becomes 

r(^(3770) -> LH) = 467+£|pCeV(±50%), (13) 

Br(V>(3770) LH) = (2.0;? j°)%(±50%). (14) 

For a small mass m c = 1.4 GeV, the LH decay width and 
branching ratio of -0(3770) can reach 805 KeV(±50%) 
and 3.5% (±50%) respectively, and this could be viewed 
as the maximum value for the LH decay of 0(3770) in 
our estimation based on the calculation at leading order 
in v 2 and next-to leading order in a s in NRQCD. 

Together with the partial decay width of 350-400 KeV 
observed for hadronic transitions and El transitions of 
the ■0(3770), the predicted annihilation (LH) decay width 
in Eq.(13) will make the total non-DD decay width of 
V>(3770) to be about 820-870 KeV for m c = 1.5 GeV, 
and 1.15-1.20 MeV for m c = 1.4 GeV. The latter may be 
viewed as the maximum value obtained in our approach 
for the total non-DD decay width, corresponding to a 
branching ratio of about 5% of the -0(3770) decay. 

In summary, we have given a rigorous theoretical pre- 
diction for the LH decay of 0(3770), based on NRQCD 
factorization at NLO in a s and LO in v 2 . By intro- 
ducing the color-octet contributions, the results are free 
from infrared divergences. We find that for the -0(3770) 
the D-wave contribution is dominant, and the effect of S- 
D mixing is very small. Numerically, our results do not 
favor either of the two experimental results measured by 
BES and CLEO collaborations. We hope our theoreti- 
cal result can serve as a clue to clarify the long-standing 
puzzle of the -0(3770) non-DD decay. We urge doing 
more precise measurements on both inclusive and exclu- 
sive non-DD decays of 0(3770) in the future. If their 
total branching ratio can be as large as 10%, it will be a 
real challenge to our current understanding of QCD, and 
new decay mechanisms have to be considered. 
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